Direct Assessment of Vorticity Alignment with Local and Nonlocal Strain Rates in 

Turbulent Flows 
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A direct Biot-Savart integration is used to decompose the strain rate into its local and nonlocal 
constituents, allowing the vorticity alignment with the local and nonlocal strain rate eigenvectors 
to be investigated. These strain rate tensor constituents are evaluated in a turbulent flow using 
data from highly-resolved direct numerical simulations. While the vorticity aligns preferentially 
with the intermediate eigenvector of the combined strain rate, as has been observed previously, the 
present results for the first time clearly show that the vorticity aligns with the most extensional 
eigenvector of the nonlocal strain rate. This in turn reveals a significant linear contribution to the 
vortex stretching dynamics in turbulent flows. 

PACS numbers: 47.27.-i,47.32.C-,47.27.De 



The alignment of vorticity with the strain rate eigen- 
vectors in turbulent flows has been a subject of consider- 
able interest over the past two decades. Since the initial 
finding [l| that the vorticity shows a preferred alignment 
with the intermediate eigenvector of the strain rate ten- 
sor, there have been numerous studies seeking to under- 
stand the reasons for this result, and various theoretical 
approaches have been proposed to explain the failure of 
the vorticity to align with the most extensional strain 
rate eigenvector. 

In this Letter, we help resolve this issue by show- 
ing that vorticity in turbulence does tend toward align- 
ment with the most extensional eigenvector of the nonlo- 
cal (background) strain, namely the strain field induced 
in the immediate region around any vortical structure 
by the surrounding vorticity outside this region. The 
anomalous alignment occurs with the eigenvectors of the 
combined strain rate, namely the sum of this nonlocal 
background strain and the local strain induced in the re- 
gion by the vorticity within it. 

Alignment of the vorticity vector u> = V x u with the 
strain rate tensor Sij in three-dimensional incompressible 
turbulent flows is ultimately responsible for the transfer 
of kinetic energy between scales, and for the nonlinearity 
in the dynamics of the underlying vorticity field. The 
inverse curl operator is the Biot-Savart integral that gives 
the velocity field u from the vorticity field uj as 



u(x) 



1 

4-7T 



w(x') X 



,/|3 



(1) 



The resulting gradients of u define the strain rate tensor 
Sij = \ (dui/dxj +duj/dxi) which, in turn, is coupled 
to the dynamics of the vorticity as 



(2) 



On the right side of ©, the magnitude of the stretch- 
ing term \SijUj\ = uj[sf (ei ■ e u ) J 1 / 2 depends on the 
strain rate eigenvalues s, and the vorticity magnitude 

1/2 

u) = (uiiUJi) , and on the alignment cosines (e* • e u ) be- 
tween the vorticity unit vector e w and the strain rate 
eigenvectors e,-. 

The strain rate eigenvalues s, can be ordered as s% > 
3 2 > S3, so that incompressibility (s% + s% + S3 = 0) re- 
quires si > and S3 < 0. The positivity of si and 
the negativity of S3 correspond, respectively, to exten- 
sional and compressional straining along the ei and e3 
directions. While the intermediate eigenvalue S2 is on 
average weakly positive in turbulent flows, the instan- 
taneous S2 can take on large positive or negative values 
[ESUl bounded only by the S\ and S3 values. The align- 
ment between the vorticity and the strain rate eigenvec- 
tors similarly determines the production rate LOiSijUJj = 
w 2 s i (e i -e w ) 2 for the enstrophy |(u;-u>). The three 
alignment cosines (e^ ■ e w ) thus play an essential role in 
the structure and dynamics of turbulent flows. 

Despite its importance, the mechanism by which the 
vorticity aligns with the strain rate eigenvectors ej is still 
not well understood. In particular, the maximality and 
positivity of si might suggest that the vorticity in @ 
would show preferred alignment with the most exten- 
sional eigenvector ei. However, since Sij on the right 
side of @ is coupled back to Wj, the resulting nonlinear- 
ity complicates any such simple alignment. Indeed, nu- 
merous studies have shown that the vorticity in turbulent 
flows instead shows a preference for alignment with the 
intermediate strain rate eigenvector e2. 

This can be seen, for example, in Figure [TJ where 
distributions of the three alignment cosines |ej • e u | are 
shown from a recent highly-resolved, three-dimensional, 
2048 3 direct numerical simulation (DNS) [H of statisti- 
cally stationary, forced, homogeneous, isotropic turbu- 
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FIG. 1: (Color online) Distributions of alignment cosines |e< • 



lence at Taylor-scale Reynolds number i?e> = 107. The 
simulation was done using a pseudospectral method with 
a spectral resolution that exceeds the standard value by a 
factor of eight. As a result the highest wavenumber corre- 
sponds to k max r)K — 10, and the Kolmogorov lengthscale 
t\k is resolved with three grid spacings. The resulting 
alignment distributions in Fig. [1] agree with those from 
lower-resolution DNS studies as well as from laboratory 
measurements [!, H, [||. In particular, the vorticity tends 
to point away from the most compressive eigenvector e3, 
namely |e3 ■ e u \ — * 0, and there is essentially no tendency 
for any preferred alignment relative to the most exten- 
sional strain rate eigenvector ej., since P (|ei • e w |) » 1. 
However, the vorticity shows a strong tendency toward 
alignment with the intermediate strain rate eigenvector 
e 2 , namely |e 2 • e w | — ► 1. 

Previous attempts to understand these alignments 
have often focused on the vorticity-strain coupling in the 
evolution equation for the strain rate tensor, namely 
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d 2 p/ 'dxidxj is the pressure Hessian. The 
nonlinear coupling between the strain rate and the vor- 
ticity is apparent in ((2]) and ((3J, where the nonlocality 
of the vorticity and strain rate co-evolution involves the 
pressure Hessian 0, H[. While progress has been made 
in understanding the vorticity alignments in Fig. [1] via 
the restricted Euler equations (e.g. H), where nonlo- 
cal effects are neglected entirely, a complete picture that 
clearly distinguishes between local and nonlocal contri- 
butions to the vorticity dynamics has remained largely 
elusive. 

Here, we forgo the use of ([3]) in representing the strain 
rate appearing in @ , instead using an integral represen- 
tation for the strain rate derived from ([1]), and then use 
this to gain insights into the effects of nonlocality and 
nonlinearity on the vorticity alignment. As shown in Q, 
from the Biot-Savart integral in (|TJ) the strain rate can 



be exactly expressed in terms of the vorticity as 

3 f fk 

Sij(x) = — (e ik irj + ntjki) — w/(x')d 3 x' , (4) 

o7T J x / T 

where r = x — x', r = |r|, and the integral is defined 
in a principal value sense. Substituting (TJJ in @ then 
gives a direct nonlocal integro-differential equation for 
the vorticity evolution as 
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which depends only on the vorticity field itself. In (H|) and 
(O , the local and nonlocal contributions to the strain rate 
and vorticity dynamics can be understood by separating 
the integration domain into a local region of radius r < R 
centered on x, and a nonlocal region that accounts for the 
rest of the domain [l(| ■ The strain rate in ([JJ then is 
the sum 



S ij -(x)=55(x) + 5^(x) 



(6) 



of the local strain rate <S^(x) induced at x by the vortic- 
ity within R, and the nonlocal (background) strain rate 
SR (x) induced at x by all the vorticity outside R. 

The background strain field Sy(x) in the vicinity of 
any local vortical structure in the turbulence is that in- 
duced by all the other vortical structures. Thus, the 
proper physical value for R used to obtain (x) should 
exclude from ((4]) essentially all the vorticity associated 
with any local vortical structure. Prior studies (e.g. fill ]) 
have shown that the characteristic radius of intense vorti- 
cal structures in turbulence is in the range r/r/K ~ 4 — 10, 
where r\K is the Kolmogorov length scale. This is con- 
sistent with the two-point vorticity correlation from the 
present DNS of homogeneous isotropic turbulence, which 
is found to decrease to 20% of its maximum value at 
r m 12rjK- This gives a physically appropriate cutoff ra- 
dius, since beyond this the vorticity becomes essentially 
uncorrelated with itself. Thus, R = Y1t\k as used herein 
excludes essentially all the local vorticity for most struc- 
tures, and thereby allows the self-induced strain field in 
the vicinity of typical vortical structures to be separated 
from the background strain field in which the structures 
reside. 

Although § developed an operator for the background 
strain rate S^(x), here we avoid the associated infinite 
Taylor series with respect to R and instead obtain the 
local strain rate S^(x) by directly integrating Q over 
the domain R centered on x. At any point x in the 2048 3 
cubic simulation domain, a smaller cubic subdomain with 
side length 2R is taken to define the local region around 
x. The local strain rate Sy (x) is obtained by numerically 
integrating Q over this subdomain. We then determine 
the nonlocal strain rate from © as Sy(x) = S'y(x) — 
£jy(x), and examine the alignment of the vorticity u>(x) 
with each of these strain rates to understand how the 
alignment in Fig. [1] arises. 
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FIG. 2: (Color online) Typical example of (a) background 
strain rate field S'S(x) and (b) local strain rate field S^x). 



Figure [2] shows an example of the resulting decompo- 
sition of the shear strain rate field S12 (x) into its back- 
ground and local fields, 5^(x) and S'^l(x). Similar local- 
nonlocal decompositions are obtained for the other strain 
rate components, and the eigenvalues and eigenvectors 
of the resulting background and local strain rate tensor 
fields are then computed. At every point x, the alignment 
cosines • e u \ of the vorticity with the background and 
local strain rate eigenvectors, denoted ef and ef respec- 
tively, are then evaluated. 

The resulting vorticity alignment distributions 
P(|ef -e u \) and P (|ef ■ e u \) are shown, respectively, 
in Figs, and [3f . From the background strain align- 
ments in Fig. [3)3, it is apparent that the vorticity is 
preferentially aligned with the most extensional back- 
ground strain rate eigenvector ef , namely ef ■ e w | — > 1. 
There is essentially no preferred alignment of the vortic- 
ity relative to the intermediate background eigenvector 
ef, since P(|ef -e w |) w 1, while the vorticity tends 
to point preferentially away from the most compressive 
background eigenvector ef , namely |ef ■ e w | — > 0. 

The alignments in Fig. ^jp, with the background strain 
rate are precisely as would be expected when the strain 
rate evolution is decoupled from that of the vorticity, as is 
essentially the case for the background strain. From $Z$i 
with the inviscid dynamics of the vorticity satisfies 
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By definition, the background strain rate $y in ([7]) is 
independent of the vorticity at x, and thus its effect on 
the dynamics of the vorticity Wi(x) is essentially linear. 
Since sf > and sf < 0, and since sf < sf , the ef- 
fect is to cause u> to rotate toward alignment with the 
most extensional eigenvector ef of Sfy The fact that 
such alignment of the vorticity is seen in Fig.[3]a suggests 
that the quasi-linear dynamics from the first term on the 
right in ([7]) plays at least a significant role in the overall 
evolution of the vorticity field in turbulent flows. In the 
terminology of She et al. , this would be referred to as 



kinematic nonlocality, as distinguished from the dynamic 
locality that was the focus of their study. 

From the vorticity alignments in Fig. \3$> with the lo- 
cal strain rate field induced in R by the local vorticity, 
u) shows substantial and essentially equal preference for 
pointing largely perpendicular to the most extensional 
and compressional eigenvectors ef and ef of the local 
strain rate S^, namely |ef ; • e u \ — > and |ef ' ■ e u \ — > 0. 
This is consistent with the fact that much of the vortic- 
ity in turbulent flows concentrates into relatively com- 
pact line-like and sheet-like structures formed by locally 
axisymmctric and planar background strain rate fields 
[H, EH- In the former case, the axisymmetric Burg- 
ers vortex is often used as an idealized representation of 
such structures, while in the latter case the planar Burg- 
ers vortex sheet provides a similar idealized representa- 
tion. In both cases, the two-dimensional local strain field 
induced by the vortical structure has large extensional 
and compressional eigenvalues with eigenvectors that are 
necessarily perpendicular to the vorticity, due to the ge- 
ometry of the structures. The remaining eigenvalue is 
zero for perfectly two-dimensional structures, and will 
be nonzero only due to small departures from strict two- 
dimensionality of the structures. Its small magnitude 
is thus nearly always between the other two eigenval- 
ues, and will therefore be the intermediate eigenvalue. 
Its eigenvector must necessarily be perpendicular to the 
other two, and so will necessarily be closely aligned with 
the vorticity itself. 
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FIG. 3: (Color online) Distributions of vorticity alignment 
cosines for (a) background strain rate eigenvectors ef , and 
(6) local strain rate eigenvectors ef . 
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FIG. 4: Distribution of background (linear) to local (nonlin- 
ear) vortex stretching ratio. 

This is precisely the alignment seen with ef in Fig.[3t, 
where the vorticity points strongly along the direction of 
the intermediate eigenvector, namely |ef -e^l — > 1. Note 
that this 'preferred' alignment of the vorticity with the 
intermediate local strain eigenvector in Fig. [3^ is not a 
result of the nonlinear dynamics from the second term on 
the right side in ([7]), but rather is a simple geometric con- 
sequence of the largely sheet-like and line-like structures 
into which the vorticity is formed [1 21 ] . 

The alignments in Fig. [3] allow the relative contribu- 
tions from the two terms on the right side of the inviscid 
vorticity dynamics in J7J to be understood. In particular, 
Fig. 2] shows the distribution of the ratio of background 
and local vortex stretching rates, namely \S^ujj\/\S^ujj\. 
It is apparent that in much of the flow this stretching 
ratio exceeds one, meaning that the linear stretching dy- 
namics produced by the background (nonlocal) strain 
field S-j in flTJ exceeds the nonlinear stretching dynamics 
from the local strain field S-j. Thus, despite the overall 
nonlinear dynamics governing the vorticity evolution in 
@ , a substantial part of the underlying dynamics is lin- 
ear and nonlocal. 

This is consistent with the alignment in Fig. of the 
vorticity with the intermediate eigenvector of the local 
strain rate S-j, for which the associated eigenvalue s§ 
has the smallest magnitude among the three local strain 
eigenvalues, and thus the associated stretching is not nec- 
essarily large. By contrast, Fig. shows that the vor- 
ticity aligns with the most extensional eigenvector of the 
background strain rate S^, and thus is stretched by the 
largest of its three eigenvalues. As a result, even when 
\SfA is smaller than \S-j\, the background stretching may 



be larger than the local stretching. This is a conse- 
quence of the fact that in ([7]) vortex stretching by the 
local strain rate is generally not favored from the stand- 
point of geometrical alignment. It is remarkable that the 
linear stretching dynamics from this background (nonlo- 
cal) strain field is comparable to the nonlinear stretching 
dynamics from the local strain field. 

Further insights into the background and local 
dynamics may be gained by conditioning the vortex 
stretching ratio on the vorticity magnitude w, as shown 
in Fig. [51 While there is a tendency towards smaller 
vortex stretching ratios for large vorticity magnitudes 
in Fig. [5j the observed dependence is relatively weak. 
This is due to the competition between increased local 
strain rate magnitude (which favors local stretching) 
and the correspondence with nearly two-dimensional 
intense vortical structures (which favors background 
stretching) for large values of u). Unraveling the in- 
dividual contributions of these two effects, as well as 
consideration of other secondary parameters such as 
the background strain persistence or the Reynolds 
number, is an important direction for future research. 
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FIG. 5: (Color online) Distribution of background to local 
vortex stretching ratio conditioned on vorticity magnitude 
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